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Abstract 

Irreversible adsorption of spheres on flat collectors having dimension d < 2 is studied. Molecules 
are adsorbed on Sierpinskis Triangle and Carpet like fractals (1 < d < 2), and on General Cantor 
Set (d < 1). Adsorption process is modeled numerically using Random Sequential Adsorption 
(RSA) algorithm. The paper concentrates on measurement of fundamental properties of cover- 
ages, i.e. maximal random coverage ratio and density autocorrelation function, as well as RSA 
kinetics. Obtained results allow to improve phenomenological relation between maximal random 
coverage ratio and collector dimension. Moreover, simulations show that, in general, most of known 
dimensional properties of adsorbed monolayers are valid for non-integer dimensions. 
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I. INTRODUCTION 



Most of the effort in the research on irreversible adsorption is focused on flat homogeneous 
surfaces as it can be further directly exploited in material sciences. However, experiments 
using fractal-like collectors were also made 1], though preparation of such surfaces is com- 
plicated. On the other hand, there are many porous media in nature and adsorption may 
play an important role there Q, [J. For example, coral fractal-like structure helps them 
to catch plankton effectively [4j. Adsorption on fractal collectors might also be applied in 
environmental protection in designing effective water or air filters jsj. 

Hard spheres packing has a rich research history, both in physics and maths. The greatest 
interest has been shown in the maximal possible packing Q] . The most common property 
of a given packing is its density. In d dimensional Euclidean space, it is given by 

e ma M = p-<d) (i) 

where p is a number of spheres in a unit volume, and 

is a volume of d dimensional sphere having radius tq. T(x) is the Euler gamma function. 

Presented work focuses on the properties of maximal random coverages, built of spherical 
particles on flat collectors having fractal dimension of a non-integer value. Earlier, similar 
studies concentrated mostly on adsorption on lattice collectors. Their popularity was pow- 



ered by a number of analytic results obtained, e.g. 
ractal collectors concentrate on diffusion properties 
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Other investigations concerning 



111 or adsorption on rough surfaces 

n 

13] . Fractal structures formed in the adsorption processes were also observed [141 ] . 



Dimensional properties of adsorption layers were studied mostly for collectors having 
integer dimensions, with special attention paid to two-dimensional case, e.g. jlolJisl]. due 
to its potential application in chemistry and material science. Here, the only analytically 
solved case is a one-dimensional problem known also as car parking problem for which 



^mas(l) — 0.748... 19[. The analysis of hard spheres random packing for 2 < d < 6, nicely 



reviewed in 



1* 



shows that maximal random coverage ratio decreases with the growth of 



collectors dimension: 



Omax(d) - 2<J 2 (3) 



with ci = 0.202 and c 2 = 0.974 
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20] . There are also some characteristics, which do not 



depend qualitatively on dimension. For example the two-point density correlation function 
'or maximal random coverages is known to have superexponential decay for large distances 



2l| 



if 2 \(»72n>)-l 

C(r) ~ W) U(r/(^)-l) J fM '' * °°' <4> 



and logarithmic singularity when particles are in touch [lg, [1 

C(r) ~ In (J- - 1) for r -)• (2r ) + and -)• # max . (5) 

The aim of this work is to check if the mentioned characteristics are valid for adsorption on 
collectors having non-integer fractal dimensions. The maximal coverages has been obtained 
here numerically by using Random Sequential Adsorption (RSA) algorithm [15| described in 
detail in Section [Til The kinetics of RSA simulation has also been of interest to the authors 
of this work, as it affects estimation of max (d). For sufficiently long simulation time the 
coverage ratio 6(t) scales with collector's dimension [15I. flf]]: 

e max - e{t) ~ r l ' d . (6) 



Here 6 max = 9{t — > 00). Importance of the above relation is increased 



)y the fact that 



it is also valid for other molecules, more complex than simple spheres [22|, |23[ . To model 
collectors having dimension from a range of d G (1,2) the Sierpinski's Triangle and, Carpet- 
like fractals are used. Adsorption for d < 1 is studied using General Cantor Set. 



II. MODEL 



Maximal random coverages are generated using RSA algorithm, which has been success- 
fully applied to study colloidal systems. It is based on independent, repeated attempts of 
adding sphere to a covering film. The numerical procedure is carried out in the following 
steps: 

i: a virtual sphere is created with its centre position on a collector chosen randomly according 
to the uniform probability distribution; 

ii: an overlapping test is performed for previously adsorbed nearest neighbours of the virtual 
particle. The test checks if surface-to-surface distance between spheres is greater than 
zero; 



iii: if there is no overlap the virtual particle is irreversibly adsorbed and added to an existing 
covering layer. Its position does not change during further calculations; 

iiii: if there is an overlap the virtual sphere is removed and abandoned. 

Attempts are repeated iteratively. Their number is usually expressed in a dimensionless 
time units: 

to = (7) 

where N is a number of attempts and Sp denotes an area of a single sphere projection on a 
collector. Typically, Sc is a collector area. Although formally the area of any fractal having 
D < 2 is zero, Sc denotes here an area of an overlapping square and is equal to (100r) 2 , 
where r is a radius of an adsorbed sphere. Collectors are modeled as a successive iteration 
of a given fractal. For each collector type at least 20 independent numerical experiments 
have been performed. Sample coverages are presented in FigJT] Dimension of a given fractal 
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FIG. 1: Maximal random coverages on first, second and sixth iteration of example fractals. 
The side length of square containing the fractal is lOOr. 



is estimated using relation 



d 



ln(p) 
hu>) 



(8) 



where p is a number of identical smaller parts created in a single iteration and s is a scale 
change. In case of Cantor Dust, for example, after one iteration four smaller squares appear 
after one iteration (p = 4) with a side length three times smaller than in the original square 
(s = 3). Therefore, dimension of Cantor Dust is ln(4)/ln(3) ~ 1.262 

Coverage ratio of a given adsorption layer is determined by using random sampling of 
collector's points and checking whether they are covered by any, previously adsorbed sphere: 

e(t) = ^ (9) 

71/ 

where n c (t) is a number of covered points after a dimensionless time t; n denotes total 
number of sampled points. Here n = 10 6 , which provides statistical error at the level of 
0.1%. 



III. RESULTS AND DISCUSSION 



A. Kinetics of the RSA 



There are two factors affecting maximal random coverage ratio obtained during simula- 
tion. First, the collector used is not an actual fractal but only a two-dimensional approxima- 
tion and the number of adsorbed spheres clearly depends on the level of this approximation 
(see FigJT]). Second, the maximal random coverage is defined for a collector having infinite 
size and is achieved after an infinite simulation time. 

To handle the first factor, the successive iterations of fractals have been used. The number 
of adsorbed particles becomes stable at about the 5 th fractal iteration (see Figj2]). Similarly, 
the coverage calculated using Eq.([9]) does not change thereafter, too (FigOJ). Therefore, 
using the 10 th iteration has been assumed to be save and to give reliable results. 

To obtain maximal random coverage ratio after an infinite time, the RSA kinetics model 
has to be used. The analysis of numerical experiments shows that RSA for spheres in integer 
dimensions obeys the Feder's law ([6]). The results presented in FigJD confirm that for almost 
all the cases analysed in this study. Additionally, the difference between observed maximal 
number of adsorbed particles for t = 10 5 and the asymptotic value does not exceed 0.2%. 
Coverage does not grow asymptotically only in case of adsorption on the sole incoherent 
collector studied, the Cantor Dust. In this case, the maximal possible number of 256 particles 
is adsorbed after relatively short time. This is due to geometrical relations between collector 
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FIG. 2: The number of adsorbed particles depends on a fractal's iteration. Iteration zero is 
a solid triangle in case of the Sierpinski Triangle and a square in all the other cases. 
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FIG. 3: Coverage ratio after t = 10 5 as a function of fractal's iteration. Iteration zero is 
a solid triangle in case of the Sierpinski's Triangle and a square in all the other cases. 



size and spheres' diameter (see FigJSJ). After four iterations, each of the remaining squares 
has a side length of 100/3 4 r ~ 1.234r and can adsorb only one sphere. Moreover, due 
to the large distance between parts of the fractal, particle adsorbed on one square cannot 
block absorption events on its neighbouring squares. Therefore, each square contains a 
single sphere. It is noteworthy that further iterations of the fractal do not affect the above 
reasoning. Therefore, after catching 4 4 = 256 spheres no further adsorption acts are possible. 
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FIG. 4: Asymptotic character of RSA kinetics for collectors of different dimensions. The 
fits correspond to Eq.(JH]). (a) d = 1.893 (Sierpinski Carpet), (b) d = 1.771, (c) d = 1.631, 
(d) d = 1.585 (Sierpinski Triangle), (e) d = 1.465 (Vicsec fractal), (f) d = 1.262 (Cantor 

Dust). 



B. The Maximal Random Coverage Ratio 

Investigating the relationship between maximal random coverage ratio and fractal di- 
mension is the main purpose of this study. Although, there are number of studies providing 
lower and upper limits for a saturated coverage, e.g. [24j, the best direct fit so far is given 
by Eq.flB]). It perfectly matches simulation data for 2 < d < 6 [18]. However, it does not 
work well for d — 1 where 6 max (l) = 0.748 and d = where 6 m ax{Q) — 1- To obtain better 
approximation, we added next order term in the nominator of relation Q: 

d + c 2 d + c 3 d 2 

Omax[d) = ~ d (10) 
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FIG. 5: Sample coverage on the Cantor Dust collector. 



Figure [6] presents plot of 9 max (d) where values for 1 < d < 2 were obtained with (jUJ) and 
whereas values for d G {2, 3, 4, 5, 6} were taken from 18]. 
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FIG. 6: Maximal random coverage ratio in different dimensions. Dots represents obtained 
data, solid line is the (j3]) fit (ci = 0.978 and c 2 = 0.581), whereas dashed line is the (fill fit 

with Cl = 0.999, c 2 = 0.416 and c 3 = 0.086. 

The fit defined by Eq. (TiUj) matches the data for d < 2 and is characterized by correlation 
coefficient R = 0.9992. Furthermore, the parameter c\ corresponds very well to random 
saturation density for d = 0. The correction given by the third term c^d 2 = 0.086 • d 2 is 
important mostly for large d values. Therefore, simulations for higher dimensions should 
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determine whether C3 should remain or not. 



Pair correlation function 



Two-point correlation functions for fractal collectors are presented in FigJTl Again, except 
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FIG. 7: Two-point correlation function for collectors of different dimensions. Insets show 
asymptotic character for particles in close proximity to each other with a fit given by (JSJ). 
(a) d = 1.893 (Sierpinski Carpet), (b) d = 1.771, (c) d = 1.631, (d) d = 1.585 (Sierpinski 
Triangle), (e) d = 1.465 (Vicsec fractal), (f) d = 1.262 (Cantor Dust). 



for the case of Cantor Dust collector, all the plots look typical; for d < 1.7, however, the 
specific fractal structure is manifested by a number of local extremes. Therefore, the pure 
superexponential decay (Hj) cannot be observed. For small values of r, the logarithmic 
singularity (JSJ) appears just as for adsorption in higher dimensions. In case of the Cantor 
Dust collector, the two-point correlation function behaviour is strictly connected with cluster 



character of the fractal. There is also no point to investigate the singularities at kissing limit 
because close packing is restricted by incoherent structure of the fractal. 



D. Below d=l 

The preliminary RSA simulation has been also performed for d < 1 using General Cantor 
Set (GCS) as a collector. The GCS dimension depends on the size of middle part removed 
in each iteration: 

**»< 4 > = ln[2/(ll g) ] < U > 
where q is the relative length of the removed fragment. For standard Cantor Set q — 1/3 
thus dcs — ln(2)/ln(3) ~ 0.631. The diameter of the collector used during simulations was 
lOOro. The dependence of 9 max on dimension value is shown in FigJHJ The data for d > 0.5 




FIG. 8: 6 max (d) for d < 1. Dots represent data and dashed line is the ( 1101) fit with 
Ci = 0.999, C2 = 0.416 and C3 = 0.086 - the same as presented in FigJHJ 



are consistent with the fit f jTOj) within an error range. For lower d, when holes in a collector 
are bigger, the geometrical relations start to play the most important role. The situation 
is somehow similar to the Cantor Dust case; however for this case, it can be proved that 
random maximal coverage ratio depends strictly on the particle diameter to collector size 
ratio. Therefore, the possibility of getting reliable data for d < 0.5 from RSA simulation is 
doubtful. There is also no point in analysing the two-point correlation function because all 
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collectors having d < 1 are incoherent and the correlation behaviour will mainly express the 
structure of a specific collector. 

IV. SUMMARY 

The RSA on fractal collectors having d < 2 was studied. For such systems, the adsorption 
kinetics obey the Feder's law fl6]). Moreover, the singularity of the two-point correlation 
function at kissing limit is logarithmic. The decay seems to be superexponential, but the 
autocorrelations for low dimensional fractals reflect also collector structure. It was one of 
the reasons why obtaining reliable results from simulation on collectors having dimension 
significantly lower than d = 1 is not possible. The maximal random coverage ratio for the 
whole range of collector dimensions is governed by the Eq. fllOp . which is a simple extension 
to the rule ([3]) introduced and analytically supported for high dimensions. 

This work was supported by grant MNiSW/0013/H03/2010/70. 
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